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Observations and interpretation of extragalactic rotational and rovibrational H2 emission are 
reviewed. Direct observations of H2 lines do not trace bulk H2 mass, but excitation rate. As 
such, the H2 lines are unique diagnostics, if the excitation mechanism can be determined, which 
generally requires high-quality spectroscopy and suitable additional data. The diagnostic power 
of the H2 lines is illustrated by two cases studies: H2 purely rotational line emission from the 
disk of the nearby spiral galaxy NGC 891 and high resolution imaging and spectroscopy of H2 
vibrational line emision from the luminous merger NGC 6240. 



1. Introduction 

Direct observations of H2 emission from external galaxies have become standard prac- 
tice in the past decade through the revolution in ground-based near-infrared instrumen- 
tation. As a result, the near-infrared H2 rovibrational lines are now readily detectable 
throughout the local universe (e.g., Moorwood & Oliva 1988, 1990; Puxley et al. 1988, 
1990; Goldader et al. 1995, 1997; Vanzi et al. 1998). More recently, the Short Wavelength 
Spectrograph (SWS) on the Infrared Space Observatory (ISO) has for the first time al- 
lowed detection of the purely rotational H2 lines in the mid-infrared spectral regime. For 
instance, the first detection (outside the solar system) of the H2 S(0) line at 28.21 /im 
was reported by Valentijn et al. (1996) from the star forming nucleus of the nearby spiral 
galaxy NGC 6946. 

The interpretation of these data is, however, far from trivial. At typical molecular 
cloud temperatures (T ~ 20 K) the upper levels of even the lowest H2 transitions are 
essentially unpopulated, and hence H2 emission is, unlike for instance CO J = 1— »-0 
emission, not a tracer of bulk molecular gas mass. Instead, an excitation mechanism 
capable of populating these energy levels is required. Furthermore, the excitation rate 
needs to be high enough to maintain an excited level population sufficient for producing 
detectable emission. Hence the H2 line luminosities measure the rate of excitation] as 
such they provide unique and highly diagnostic information, that cannot be obtained in 
any other way. In addition, since the H2 lines are forbidden quadrupole transitions with 
very small Einstein A values, the lines are optically thin; hence there is, in contrast to 
the situation with most other molecular lines, no need to solve a complicated, geometry- 
dependent radiative transfer problem for a physical interpretation of the H2 lines. The 
near- infrared H2 lines will of course suffer from extinction by dust, but this effect can 
usually be quantified adequately by using suitable ratios of hydrogen recombination lines, 
[Fen] lines or H2 rovibrational lines in the same spectral range, with accurately known 
intrinsic flux ratios. 

The situation is complicated, however, by the fact that H2 excitation can be brought 
about by a variety of mechanisms, which may all play a role. For instance, in a starburst 
galaxy, H2 emission may be generated by UV-pumping in starforming regions or by shocks 
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due to supernova remnants or outflows; moreover, these processes are expected to occur 
together in the same small volume occupied by the starburst, and thus a combination of 
excitation mechanisms may be expected at every position. Generally, this combination 
will be difficult to separate (e.g., NGC253 observations: Forbes et al. 1993; Engelbracht 
et al. 1998). In addition, shocks due to large-scale streaming motions in spiral arms, 
bars (favoured for the barred Seyfert 2 NGC 1068 by Tacconi et al. 1994) or merger- 
driven flows (e.g., in NGC 6240, Van der Werf et al. 1993; see also §|^) may play a role. 
Furthermore, X-ray excitation may be produced by multiple supernova remnants, or, if 
present, by an active galactic nucleus (proposed for e.g., Cyg A by Ward et al. 1991 and 
for Cen A by Bryant & Hunstead 1999) or a cooling flow (Jaffe & Bremer 1997). In the 
absence of sufficient spatial resolution for separating the various excitation mechanisms, 
spectral diagnostics must be used. Models for UV-pumped, shocked or X-ray excited 
H2 emission have reached considerable predictive power; however, the densities of the 
emitting regions are often sufficiently high to thermalize the relevant level populations, 
quenching the typical signatures of non-thermal excitation processes. The only accessible 
diagnostics are then fluxes of lines with very high critical densities. However, these lines 
are intrinsically faint, and accurate fluxes for such lines require long integrations with 
large telescopes. The complexity of the problem is well illustrated by the H 2 emission of 
the merging system NGC 6240, which has been attributed to X-ray excitation (Draine & 
Woods 1990), UV-pumping (Tanaka et al. 1991), shocks (Van der Werf et al. 1993) and 
formation pumping (Mouri & Taniguchi 1995). However, all of the non-thermal excitation 
processes relied on poorly measured fluxes of faint lines. More recent spectroscopy of 
NGC 6240 by Sugai et al. (1997) showed that only shock excitation can account for the 
H2 rovibrational spectrum (see also § ||) . 

Generally a combination of accurate multi-line spectroscopy and suitable additional 
information such as high resolution spatial information is needed for a proper analysis of 
the dominant excitation mechanism and hence a physical analysis of the H2 emission. A 
complicating factor when combining H2 rovibrational and purely rotational lines is the 
fact that the dominant excitation mechanisms of these lines may be different, because 
of the different excitation requirements of these transitions: while the lowest rotational 
lines require T > 100 K and are excited at moderate densities (nn 2 > 10 2 cm -3 ), the 
rovibrational lines require T > 2000 K and nn 2 > 10 4 cm" 3 . The following sections will 
therefore present two case studies. First an analysis of the extended purely rotational H2 
emission from the nearby edge-on spiral galaxy NGC 891 is given (§§). Then, the more 
extreme conditions in the merger NGC 6240 are discussed, using new high resolution 
near-infrared data. The general implications of these two cases are discussed in §||. 

2. Extended H2 rotational line emission in the disk of NGC 891 

Valentijn & Van der Werf (1999a) have observed seven positions in the disk of the 
nearby (distance D = 9.5 Mpc) edge-on spiral galaxy NGC 891 with the ISO SWS. The 
positions observed include the nucleus and positions spaced along the disk to galactocen- 
tric distances R of 8 kpc south of the nucleus and 11 kpc north of the nucleus. At the 
11 kpc north position, the 12 CO J = 1^0 (Garcia-Burillo et al. 1992; Scoville et al. 1993; 
Sakamoto et al. 1997), and dust emission (Israel et al. 1999) are barely detected. With 
ISO, the H 2 S(0) (28.21 ^im) and S(l) (17.03 /im) lines were detected at all positions. 
These are the first detections of these lines outside starburst or active nuclei, with the 
exception of the detection in the disk of NGC 6946 reported by Valentijn & Van der Werf 
(1999b). 

The simplest analysis of these data would assume that the emission is fully thermalized 
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Figure 1. H2 S(l) (upper panel) and S(0) (lower panel) lines from R — 8 kpc south in the 
disk of NGC 891, with fits for a point source and an aperture-filling source as indicated. 



(i.e., the high-density limit is assumed) and arises from an isothermal gas layer with 
an ortho-para ratio of three, in agreement with the statistical weights of the ortho- 
and para- varieties. Under these assumptions the S(1)/S(0) ratio yields the temperature 
of the emitting region, which can be combined with the observed fluxes to give the 
column density of emitting H 2 , averaged over the SWS aperture. However, these simple 
assumptions lead to unacceptable results. For instance, at the position 8 kpc north 
(a typical "disk" position), the S(0) and S(l) data then imply 7V(H 2 ) = 2.7 • 10 22 at 
T = 76 K. This result is incompatible with the CO J = 2— >1/J = 1— >0 line ratio of 
0.75, at this position (Garcia-Burillo et al. 1992), which strongly rules out a large mass 
of dense molecular gas at T — 76 K. This argument shows that the simple assumptions 
used above do not suffice. 

The discrepancy can be removed by allowing lower densities, and therefore subthermal 
excitation of the H2 J=3 level, giving rise to a higher implied temperature for a given 
S(1)/S(0) ratio. Because the temperatures involved are still much lower than the upper 
levels of the transitions involved, the implied H2 mass will be extremely sensitive to 
the adopted temperature, and thus strongly decrease as lower densities are allowed. 
An additional advantage of this procedure is that the CO J— 2 level will now also be 
subthermally populated, so that CO J = 2— >1/ J = 1— >0 line ratios lower than the 
thermal value can be tolerated in the analysis. Indeed the multi-level multi- isotope CO 
data by Garcia-Burillo et al. (1992) indicate dominant densities of riH 2 ~ 10 3 cm~ 3 in 
the nuclear region and riH 2 ~ 200 cm~ 3 throughout the disk of NGC 891. Adopting these 
densities, the ISO data imply H2 column densities varying from 10 to 30% of the value 
derived from CO, at temperatures of 120 to 130 K, throughout the inner disk (R < 8 kpc) 
of NGC 891. 
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At the outer disk positions (R = 8 — 11 kpc) however, this solution is not adequate, 
since the S(0) and S(l) line profiles at these positions are significantly different. The 
SWS uses an aperture much larger than the diffraction beam and hence the observed line 
width depends on the extent of the emission region; the S(0) line shows the broad profile 
characteristic of aperture-filling emission, while the S(l) profile is narrow, indicating 
emission from a region much smaller than the aperture (Fig. [I]). Thus the S(0) and S(l) 
lines in the outher disk of NGC 891 must originate in separate regions: a warm component 
(T > 130 K) dominating the S(l) emission and located in isolated regions in the disk and 
a separate more pervasive component, dominating the S(0) line. The latter component 
may contain very significant mass, depending on its temperature, and the implications of 
this possibility for the mass budget at the outer positions have been explored by Valentijn 
& Van der Werf (1999a). However, a solution where the component dominating the S(0) 
line is cool (T ~ 90 K) is problematic, since the heating required for maintaining a very 
large mass of H2 at T ~ 90 K cannot be provided. In thermal equilibrium the heating 
rate should equal the cooling rate, which is dominated by [C 11] 158 /Ltm emission (which 
at T ~ 90 K is an extremely efficient cooler, since the upper level of the 158 /j,m transition 
is at 91 K), with possible contributions from H 2 rotational lines, CO rotational lines, and 
[Cl] emission. The [C11] 158 /im emission from NGC 891 has been observed by Madden 
et al. (1994), and it is easily verified that in the outer disk the [Cli] emission dominates 
the H2 emission (SWS data discussed here), CO emission (barely detected in the outer 
disk) and [Cl] emission (Germ & Phillips 1997). If all available carbon (assumed to have 
solar abundance) is in the form of C + , the [C 11] emission allows a column density of at 
most N(R 2 ) ~ 7 • 10 21 cm" 2 at T - 90 K in the outer disk, averaged over the SWS 
aperture, while iV(H2) ~ 10 23 cm -2 would be required to produce the observed S(0) line 
flux. 

This problem can be solved by relaxing the final assumption: that of an ortho-para of 
three. A lower ortho-para ratio raises the implied temperature and lowers the implied 
mass in the same way as a lower density does (as discussed above). Assuming an ortho- 
para ratio of unity, the implied mass is lowered sufficiently that the required heating can 
be accounted for. 

This analysis thus favours an interpretation in which the H2 rotational emission arises 
in low-density (nn 2 ~ 200 cm -3 ), warm (T ~ 120 K) molecular gas with an ortho- 
para ratio of about unity, which pervades the disk of NGC 891 at least to the end of 
the detectable CO disk and dominates the S(0) emission; it contains 10 to 30% of the 
mass implied by CO observations (lower fractions will result if an ortho-para ratio below 
unity is adopted); concentrations of warmer gas (plausibly identified with active star 
forming regions) are ubiquitous in the inner disk and rarer in the outer disk and give 
rise to the S(l) emission. Given the density and temperature of the extended warm 
gas, this component can most likely be identified with extended low-density photon- 
dominated regions (PDRs) which form the warm envelopes of giant molecular clouds, 
heated by the local interstellar radiation field. Observations of Galactic molecular cloud 
edges in the 21 cm Hi line reveal the presence of such warm envelopes by a "limb 
brightening" in Hi (Wannier et al. 1983; Van der Werf et al. 1988, 1989). Modeling of 
molecular cloud envelopes by Andersson & Wannier (1993) shows that the temperatures 
and densities estimated here are typical for such regions, especially in the zone where 
the transition from molecular to atomic hydrogen takes place. This interpretation of the 
present H2 data in terms of extended diffuse PDRs in cloud envelopes is corroborated by 
the excellent numerical agreement with column density estimates based on [C 11] 158 pm, 
the principal coolant for such regions. In addition, the fact that the implied ortho-para 
ratio is significantly lower than three also points to a low-density PDR origin of this 
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emission, since shocks or hot, high-density PDRs would produce an ortho-para ratio of 
three (the high temperature thermal value) by spin exchange reactions with H and H + . 
In summary, the H2 rotational lines in the disk of NGC 891 arise in warm, extended 
molecular cloud envelopes; these envelopes provide a physical link between the giant 
molecular clouds and the atomic medium in which they are embedded (Chromey et al. 
1989). The H 2 rotational lines provide unique diagnostics of these regions. 

3. Vibrational H2 emission in the nuclear region of the luminous 
merger NGC 6240 

The vibrational H2 emission of NGC 6240 has attracted attention because of its ex- 
traordinary luminosity: 7 • 10 7 L Q in the H2 v — 1— >0 S(l) line alone (for Ho = 
75 km s _1 Mpc -1 and with no correction for extinction). This line contains 0.012% 
of the bolometric luminosity of NGC 6240, which is considerably higher than any other 
galaxy (Van der Werf et al. 1993). Together the vibrational lines may account for 0.1% 
of the total bolometric luminosity. 

Imaging of the H2 v — 1— »0 S(l) emission from NGC 6240 has shown that the H 2 
emission peaks between the two remnant nuclei of the merging system (Van der Werf 
et al. 1993). This morphology provides a unique constraint on the excititation mech- 
anism, since it argues against any scenario where the excitation is dominated by the 
stellar component (e.g., UV-pumping, excitation by shocks or X-rays from supernova 
remnants). Instead, the favoured excitation mechanism is slow shocks in the nuclear gas 
component, which, as shown by recent high resolution interferometry in the CO J = 2— >1 
line (Tacconi et al. 1999), also peaks between the nuclei of NGC 6240. A multi-line H 2 
vibrational spectrum (Fig. |^) confirms this excitation mechanism (Van der Werf et al. 
2000), in agreement with more limited spectroscopy by Sugai et al. (1997). The inter- 
pretation in terms of slow shocks also naturally accounts for the high ratio of H2 line 
emission to infrared continuum emission, which is a characteristic of such shocks (e.g., 
Drainc et al. 1983). 

What is the role of these shocks? In the shocks mechanical energy is dissipated and 
radiated away, mostly in spectral lines (principally H 2 , CO, H 2 and [Oi] lines). This 
energy is radiated away at the expense of the orbital energy of the molecular clouds in 
the central potential well. Consequently, the dissipation of mechanical energy by the 
shocks will give rise to an infall of molecular gas to the centre of the potential well. 
Therefore, the H 2 vibrational lines measure the rate of infall of molecular gas into the 
central potential well in NGC 6240. This conclusion can be quantified by writing 

-^rad = -^disj (3-1) 

where L r ad is the total luminosity radiated by the shocks and Ldis the dissipation rate 
of mechanical energy, giving rise to a molecular gas infall rate Afn 2 given by 

L dis = ^Mn 2 v 2 , (3.2) 

where v is the circular orbital velocity at the position where the shock occurs. 

Using a i^-band extinction of 0™15 (Van der Werf et al. 1993), the total luminosity 
of H2 vibrational lines from NGC 6240 becomes 7.2 • 10 8 L©; inclusion of the purely 
rotational lines observed with the ISO SWS approximately doubles this number, so that 
L rad = 1.5-10 9 L Q . 

In order to use this number to estimate -Mh 2 , h is necessary to establish more accu- 
rately the fraction of the H2 emission that is due to infalling gas. Observations with 
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Figure 2. H+K-b&nd spectrum of NGC6240 showing H 2 lines up to the v = 1— >0 S(9) line 

(Van der Werf et al. 2000). 



NICMOS on the Hubble Space Telescope (HST) by Van der Werf et al. (2000) provide 
the required information (Fig. ||). The NICMOS image shows that the emission consists 
of a number of tails (presumably related to the superwind also observed in Ha emission) , 
and concentrations assocated with the two nuclei, and a further concentration approx- 
imately (but not precisely) between the two nuclei. The relative brightness of the H2 
emission from the southern nucleus is deceptive, since this nucleus is much better centred 
in the filter that was used for these observations than the other emission components, 
in particular the northern nucleus. Taking this effect into account, it is found that 32% 
of the total H2 flux is associated with the southern nucleus, 16% is associated with the 
northern nucleus, and 12% with the component between the two nuclei, the remaining 
40% being associated with extended emission. Using inclination-corrected circular ve- 
locities (from Tecza 1999) of 270 and 360 km s _1 for the southern en northern nucleus 
respectively, and of 280 km s~ x for the ce ntra l component (Tacconi et al. 1999), the mass 
infall rates derived using Eqs. ( [jj| ) and (3.2) are 80 M© yr _1 for the southern nucleus, 
22 Mq yr _1 for the northern nucleus and 28 M yr _1 for the central component. 

The derived molecular gas inflow rate to the two nuclei is remarkably close to the 
mass consumption rate by star formation of approximately 60 Mq yr , indicating that 
the H2 emission from the nuclei directly measures the fueling of the starbursts in these 
regions. The gas inflow towards the central component is more remarkable, since this 
component is not associated with a prominent stellar component. The gravitational 
potential at this position is therefore most likely dominated by the gas itself, which is 
also indicated by high resolution interferometry in the CO J = 2— *1 line (Tacconi et al. 
1999). The absence of prominent star formation at this position shows that the central 
gas component is gravitationally stabilized, possibly by a high local shear. However, as 
the central gas column density increases, the shear must eventually be overcome and 
given the high gas density (and hence short free-fall time) an explosive starburst will 
result. In that stage NGC6240 will become a true ultraluminous infrared galaxy. 
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Figure 3. High resolution imaging of NGC 6240 in the H 2 v = 1— >0 S(l) line and the 2.15 um 
continuum with NICMOS/HST (Van der Werf et al. 2000). 



4. General implications 

The two case studies discussed above illustrate the unique diagnostic power of H 2 
emission lines provided that the excitation mechanism can be determined. The extended 
PDR emission detected in NGC 891 reveals the presence of a widespread fluorescent 
component, which is probably a common feature of star forming galaxies. For instance, 
Pak et al. (1996) have detected extended diffuse H2 vibrational emission from the inner 
Milky Way and argue that this emission is UV-excited; similarly, Harrison et al. (1998) 
argue for purely fluorescent extended H2 emission in the moderate luminosity starburst 
galaxy NGC 253. 

The configuration of shocked H2 emission from a pronounced molecular gas concentra- 
tion between the two nuclei as found in NGC 6240 is not unique either: it is also found in 
Arp220 (Van der Werf & Israel 2000), and several other mergers. Perhaps the best case 
in point is the well-known merging "Antennae" system (NGC 4038 — 4039) where pro- 
nounced CO emission is found in the interaction zone between the two nuclei (Stanford 
et al. 1990), which is the site of vigorous obscured star formation (Mirabel et al. 1998). 
These geometries may be due to the fact that the gas components are dissipative, and 
thus merge on a shorter time scale than the dissipationless stellar components, which 
merge by the slower process of dynamical friction. The H 2 emission thus provides unique 
insight into the physics of these gas-rich mergers. 

I would like to thank my collaborators Frank Israel, Alan Moorwood, Tino Oliva, 
and Edwin Valcntijn for discussions on this subject, and Guido Kosters for reducing the 
NICMOS data of NGC 6240. 
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